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SUMMARY 

A constitutive lubricant friction model for rolling/sliding concentrated contacts 
such as gears and cams was developed, based upon the Johnson and Tevaarwerk fluid 
rheology model developed earlier. The friction model reported herein differs from 
the earlier rheological models in that very large slide to roll ratios can now be 
accommodated by modifying the thermal response of the model. Also the elastic 
response of the fluid has been omitted from the model, thereby making it much simpler 
for use in the high slide to roll contacts. The effects of this simplification are 
very minimal on the outcome of the predicted friction losses (less than 1%). In 
essence then the lubricant friction model developed for the high slide to roll ratios 
treats the fluid in the concentrated contact as consisting of a non-linear viscous 
element that is pressure, temperature and strain rate dependent in its shear 
response. The fluid rheological constants required for the prediction of the 
friction losses at different contact conditions were obtained by traction 
measurements on several of the currently used gear lubricants. In these tests the 
contact pressure was varied from 1 to 1.9 GPa, the rolling speed from 10 to 30 m/sec 
and the fluid inlet temperature from 30 to 80°C. The range of these test conditions 
corresponds to those encountered in gears and cams. The fluid constants were 
extracted from the measured traction results by fitting the rheological model against 
it. 

An example calculation, using this model and the fluid parameters obtained from 
the experiments, shows that it correctly predicts trends and magnitude of gear mesh 
losses measured elsewhere for the same fluids tested here. 
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NOMENCLEATURE. 

Below follows a list of symbols used in the text and their units. 


Sym. 

DESCRIPTION 

Units 

a,b 

Semi Hertzian contact size in the x and y direction 

[m] 

A 

Fluid viscosity temperature parameter 

[°C] 

B 

Fluid viscosity pressure parameter 

[°C/Pa] 

^s 

Specific heat of the disc material 

[J/kg.°Cj 

C 

Shear stress temperature constant 

[°C/Pa] 

Cl, 2 

Number of counter pulses 

H 

^d 

Operating centre distance of gear set 

Em]’ 

d p 

Pressure solidification temperature 

[°C] 

Dy 

Viscosity solidification temperature 

[°C] 

E>s 

Shear stress solidification temperature 

[°C3 

De 

Deborah number 

[-] 

E 

Viscosity constant for non linear thermal model 

[Pa.sec] 

E f 

Composite elastic modulus for the disc material 

[Pa] 

EO 

Experimental traction viscosity constant 

[Pa.sec] 

El 

Experimental traction stress constant 

[°C/Pa] 

E2 

Experimental time delay constant 

E-] 

F( ) 

Dissipative function for traction model 

[sec" 1 ] 

Fs 

Dimensionless thermal resistance of the film 

C-] 

F x 

Contact force in the x direction 

EN] 

F y 

Contact force in the y direction 

[N] 

F z 

Normal force on the contact 

EN] 

Fi 

Contact force per unit length 

EN/m] 

G 

Fluid shear modulus (uncorrected) 

[Pa] 

Ge 

Johnsons elasticity parameter 

[-] 

0 

g 

Fluid shear strain rate 

[1/sec] 

h 

Central film thickness in contact 

Em] 

I 

Gear ratio 

H 

J1 

Dimensionless longitudinal slip variable 

[-] 

J4 

Dimensionless longitudinal traction variable 

H 

J4 e 

Elastic stress portion of J4 

[-] 

J4 p 

Plastic stress portion of J4 

[-] 

J4 t 

Thermal dimensionless longitudinal traction 

[-1 

k 

Contact aspect ratio (b/a) 

[-] 

k f 

Thermal conductivity of fluid 

[N/sec °C] 

k s 

Thermal conductivity of disc material 

[N/sec °C] 

k' 

Fluid thermal conductivity constant 

[N/sec.°C] 

k" 

Fluid thermal conductivity constant 

[°C] 

K 

Calibration constant in side slip measurement 

[m] 

K t 

Thermal resistance of gears to ambient 

[°C.sec/N.m] 

N P 

Number of pinion teeth 

[-] 
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Ng Number of gear teeth [-3 

p pressure [Pa] 

P Mean Hertz contact pressure [Pa] 

P 0 Hertzian peak contact pressure [Pa] 

P e Psuedo Peclet number ,H 

P r Reduced Hertz pressure [Pa] 

q Thermal heat flux due to traction [N/m sec] 

Q Total heat flow due to traction [N.m/sec] 

R xe Equivalent radius of curvature in x direction [m] 

R e Equivalent radius of curvature for discs [mj 

Rbp Base radius of the pinion [m] 

Rbg Base radius of the gear [m] 

S Auxiliary variable used in elastic/plastic model H 

S Roller slip H 

T Fluid shear stress [P&3 

T s Non linear shear stress parameter [Pa] 

T c Limiting shear strength of fluid [Pa] 

T 0 Limiting shear strength at inlet temperature [Pa] 

t time [sec] 

U Rolling speed [m/sec] 

AU Slip velocity in x direction [m/sec] 

V Viscosity [Pa.sec] 

VO Viscosity constant [Pa.sec] 

Av Side slip velocity of the discs [m/sec] 

Wp Angular velocity of the pinion [rad/sec] 

X Dimensionless mesh point distance [-3 

Ax Small displacement of the displacement transducer [m] 

YI Inlet shear heating factor H 

Z$ Dimensionless inlet temperature H 

GREEK SYMBOLS. 

0 Side slip angle [rad] 

9 Temperature of the fluid [°C] 

0 C Shearplane temperature [°C] 

0 O Inlet or blank temperature [°C] 

0j Inlet temperature of oil into gear box [°C] 

© Operating pressure angle for gearset [rad] 

$ Time delay parameter H 

p s Density of the disc material [kg/m"3] 

X t Hertzian contact shape parameters [-] 
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!-0 INTRODUCTION ' they ape one of the most common 

Gears play a major role in todays soc y ^ important that they 

forms used for the transmission o pow . power that is lost in 

operate under conditions of low power loss and 1 ong M* ^ ^ ^ due t0 the 

a gear system may to the momentum transfer of oil that 

churning of the gears in the oil, the rotation 0 f the gears in the oil air 

is flung off the gears, the losses ue between the gears, and losses due to 

mixture, the losses at the contact °c ^ ^ concentrate on one particular 
support bearings and seals, b th P To achieve this we wiU use the 

aspect of the losses name ye oreviously, but modified such that the 

— -r,“r;,:r= * « -«• <• - “** *- a 

concentrated contacts. 

1-1 MESH LOSSES IN GEARS the contac t between two truncated 

s „„ r„ o, ..c 

sections ot tomtom ™ mov> , met,, to trom 

equivalent roller would vary as contact. The power transmitted 

the point of initial contact to the m load betwee n the two teeth in contact, 

from one gear to the other results surfaces and the generation of a high 

This causes a local elastic deforma ion ° introduced into the mesh to 

interface stress. For most gears a ^‘provide for some form of 

prevent surface damage to reduce th ftic , ^ ^ int<j ^ contact zone an d 

cooling. The rolling motion of the tee contact area. One of the mesh 

a thin layer of this fluid will separa e thought of as the power 

losses results from this film forming action andca^ ^ from fiPst 

required to pump this thin layer of th g proportional to the thickness 

principles it can be shown ‘^f'^^^t oftL caUed the rolling friction 
of the fluid layer in the contact. 

loss. * of each tooth is dictated by its own 

At the point of contact the surface ve For aU but the pitch 

base geometry and the distance away from he pitch p in the 

point these velocities are different so that 1S encountered due to this shear 
contact is sheared. The level o resis ® n ioss< The magnitude of this level of 
gives rise to the contact sliding r • sliding velocities, the local 

friction depends on the amount of ? he nat ure of the 

temperatures, the roUing velocities, behavior of the fluid in the 

lubricating fluid. The actual shear or r F ^ ^ researc h in recent times, and is 
contact gap has been the sub]eCt g xtreme conditions in the contact zone. Close 
very complex in nature because through the contact gap reveals that it 

examination of the fluid history atmospheric to possibly several Giga Pascal 

experiences a sudden pmsi^ V ' se ' of ^ fluid in the contact results 

in a time period of 1 to .1 m/sec. 
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in a resisting shear stress through the thickness of the fluid layer. This in turn 
will lead to heat generation and from simple calculations, temperatures in the center 
of the film can easily reach several hundred degrees Centigrade. 

To be able to predict the gear mesh losses it is therefore imperative that we fully 
understand and model the behaviour of the fluid under the conditions as present in 
the gear mesh. Yet to study the rheological properties of the fluid under these 
conditions precludes the use of most of the conventional instruments used for steady 
state measurements. In fact the only suitable type of instrument for the study is a 
disc traction or friction test machine where most of the conditions are the same or 
similar to those in gear teeth in contact. 

1-2 PRIOR INVESTIGATIONS. 

The topic of power losses in gears has concerned researchers from the very 
beginning of the introduction of gears. Notable contributions are those by 
Buckingham [1963], Merrit [1972], Shipley [1962], and Chiu [1975]. A very complete 
review of the of friction prediction techniques available to gearing engineers is 
given by Martin [1978]. Anderson and Loewenthal [1981] did a very detailed review of 
the methods used by the above references and included their own. They state that the 
major differences between the various methods is in the treatment of the sliding 
losses in the mesh. 

Anderson and Loewenthal [1981] tried to use the various gear efficiency models on 
data from the literature. This data is however quite scarce and also it was not 
possible to use the most recent traction theories for the sliding loss because no 
traction data was available for the fluids used. The comparison between the 
predicted and measured results from the Anderson and Loewenthal [1981] investigation 
is therefore subject to the possible errors introduced due to the lack of good 
traction data for the lubricants. The opportunity to remove some of the uncertainty 
from the prediction methods came with the work of Mitchell and Coy [1982] , where the 
effect of the various lubricants on the power efficiency of a helicopter transmission 
were measured. The data was taken under very carefully controlled conditions and 
with well known fluids. Upon the request of Anderson and Loewenthal the transmission 
efficiency was measured for a traction fluid. Traction fluid is normally used in 
traction drives because it has much higher slip traction characteristics than 
conventional fluids. This fluid indicated a very much lower power efficiency than 
other fluids of comparable viscosity. It should also be mentioned here that the 
observed trends in efficiency with temperature were unexpected in that some of the 
fluids showed a higher loss at higher fluid temperatures. This indicates that the 
losses cannot simply be viscosity related. 

An attempt by Coy, Mitchell and Hamrock [1984] to explain the trends in the measured 
efficiencies on the basis of viscosity and pressure viscosity data was not 
successful. They could not predict the temperature efficiency trend or the correct 
order for the efficiency as measured for the fluids. It was clear that some other 
parameters played a role in the gear efficiency. 
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1-3 CURRENT RESEARCH PROGRAM. 

As pointed out by Anderson and Loewenthal [1981] the traction characteristics of 
the lubricating fluid are very important in predicting overall mesh efficiencies. 

Since the fluids used in the Mitchell and Coy [1982] were all readily available, it 
was decided that these fluids would be tested for their traction characteristics and 
that the resulting experimental traction data would be analyzed on the basis of 
recent thermal traction theories. In order to shed some light on the observed 
efficiency trends as reported by Mitchell and Coy [1982] it was agreed that a select 
number of the fluids as used there would be tested for their traction. The selected 
fluids are; 

Designation in Commercial name 

Mitchell & Coy 


F 

K 

E 

A 

I 

X 

Y 


Mobil RL-71 
Roy co 555 
Syntech NS-6774 
Kendall ATF-220 
Mobil Jet II 

Santotrac 50 (not shown in Mitchell and Coy) 
Sunoco TDF88 (not shown in Mitchell and Coy) 


The efficiencies as measured by Mitchell and Coy [1982] are shown in Fig 1 below as a 
function of inlet oil temperature. Note the trends of the efficiency with inlet 
temperature for the fluids E and X. 

The investigation reported herein was performed by Transmission Research, Inc of 
Cleveland, Ohio, under contract to the NASA Lewis Research Center. The NASA 
technical project manager was D.A. Rohn of the Advanced Concepts and Mechanisms 
Section. 
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2-0 EXPERIMENTS 

The various traction experiments were carried out on an existing twin disc test 
facility, shown in figure 2-1, 2-2 and 2-3. The test facility was also used for the 
traction test results as reported in Tevaarwerk [1985]. This test facility was 
modified such that it would be capable of measuring the longitudinal slip occurring 
in the contact, and to extend the range of sideslip. Traction curves are obtained by 
using a combination of longitudinal and the side slip techniques; so that large 
traction transfer can be measured without the need for a large motor generator set. 

It also has the advantage that the effects of bearing friction and rolling traction 
can be corrected for by the simple measurement of the longitudinal slip. This is 
particularly important when the expected levels of fluid friction are low, as is the 
case with some of the lubricant tested here. With longitudinal traction measurements 
the correction for bearing friction and rolling traction is often of the same 
magnitude as the measured friction with these lubricants. 

2-1 DESCRIPTION OF TWIN DISC MACHINE. 

For an extensive description of twin disc traction testers the reader is referred 
to the literature; Smith [1965] and Johnson and Roberts [1974]. Basically the 
machine consists of two discs, called the upper and lower disc. The lower disc is 
mounted in rolling element bearings through shafts and the only degree of freedom is 
one of rotation about this axis. The lower disc always has a transverse radius of 
curvature of infinity. The upper disc is contained in bearings that are mounted in 
the upper assembly. This upper assembly is suspended with elastic hinges such that 
only direct normal motion or axial motion is possible. The assembly will however 
always stay horizontal. The upper disc (or toroid) has curvatures so that the 
desired contact geometry is arrived at. The upper assembly is constructed such that 
the toroidal axis can be tilted relative to the horizontal plane so as to introduce 
spin on the contact. It can also be skewed about the normal to the contact to 
introduce a side slip velocity. To introduce longitudinal slip it is possible to 
attach a electro mechanical brake to the toroid and extract torque from this system. 
Control of the amount of slip resulting from this however is difficult because of the 
nature of the traction slip response of the fluids. 

In order to achieve the range of pressures required for the traction data a set of 
discs with a nominal aspect ratio of 1.7 is used. These discs are made of AISI-01 
steel, hardened to 7.00 GPa, ground and polished to a surface finish of less than .05 
jim RMS and with an out of roundness error of less than 5 pm. Between tests the 
discs are inspected for surface damage and if needed, reground and polished to bring 
them back up to specifications. The required normal load, obtained by a dead loading 
technique, can be calculated from the Hertz theory for elastic bodies in contact. 

The maximum contact normal stress P 0 is given by; 

, 1 3 / 3F z E' 

0 2 wpX y E^ S 


(2-1) 



Page 5 


where F z = contact normal load [Nj 

P 0 = Hertzian contact stress [Pa] 

E' ,= composite elastic modulus [Pa] ; ( 231 GPa for steel) 
equivalent disc radius [m3 

= (1/Rx + 1/Ry)"" 1 

V , X = Hertzian contact shape factors [-3 

2-2 INSTRUMENTATION OF THE TRACTION TESTER. 

By suitably instrumenting the disc machine the relevant experimental parameters can 
be measured. In the experiments reported here the sideslip, sideways traction force, 
toroid surface temperature, rolling velocity, the amount of longitudinal slip, and 
the degree of asperity contact are measured. The technique of measuring each of 
these variables will be discussed next. 

2-2-1 THE MEASUREMENT OF THE TRACTION FORCE. 

A ring dynamometer type load cell is used to measure the side slip force of the upper 
toroid assembly. The electrical signal from the load cell is conditioned for noise 
and amplified using common mode rejection techniques. The gain on the amplifier is 
adjusted so that a good range on the signal is measured for each test. Calibration 
of the load cell is done in situ by dead loading. This calibration is checked 
periodically. 

2-2-2 THE MEASUREMENT OF SIDE SLIP. 

The skew angle is measured by using a direct current displacement transducer on the 
upper assembly and thereby measuring the rotation angle of this assembly. This skew 
angle gives the amount of side slip/roll ratio through the relationship; 

(2-2) Av/U= tan (j^) 

where p = side slip angle [rad] 

4V = side slip velocity [m/sec] 

U = rolling velocity [m/sec] 

By measuring the amount of skew with the displacement transducer the side 
slip/roll ratio is obtained directly through; 

(2-3) Av/U=K^x 

where ^x is the displacement of the transducer and K a scale factor. The electrical 
output of the displacement transducer is filtered in an R-C network to provide a 
low-pass signal. The maximum frequency response of the R-C network is .1 sec. The 
displacement transducer is calibrated by rotating the upper assembly through a known 
angle and then calculating the amount of side slip for this angle. 
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2-2-3 THE MEASUREMENT OF ROLLING VELOCITY AND LONGITUDINAL SLIP. 

The rotational velocity of both the disc on the toroid are measured by using a MC 
6840 frequency counter. The proximity probes for the control of the counters are 
mounted near the support bearings, see Fig. 2-3. Each shaft had a simple protrusion 
on it that would produce a pulse once per revolution. This pulse is used to turn 
countdown counters on when controlled to do so. The next successive pulse from the 
proximity probe would stop the countdown. When each of the two counters had 
completed their cycle a flag is set and the register contents are read. From the 
difference between the original and the final contents of the registers the amount of 
time for a single revolution of each shaft is calculated. From this the angular 
speed for each shaft is calculated, and hence the peripheral velocity of each roller 
can be calculated if the roller radius is known. The countdown frequency on the 
counters is selected such that an accuracy of at least 1 : 10000 is obtained. By 
employing this method the slip of the rollers could be measured every 5 or six 
revolutions of the shafts. From the count down and the undeformed roller radii the 
slip is calculated as follows; 


(2-4) 


(U2 - Ul) 
(U2 + Ul) 



S) / ( H + 

Cl R1 


C2 

Cl 


) 


[-] 


This can be approximated if the amount of slip is small by the following; 


(2-5) 


S = 


R2 C2 
R1 " Cl 


where S = roller slip [-3 

U = rolling velocity [m/sec.] 

C = nunber of count down pulses [-3 
R = roller radius [m3 

Suffix 1 denotes the toroid and suffix 2 denotes the disc. 

It should be stressed that the actual amount of slip will be somewhat different 
when two different roller materials are in contact because of the deformation of the 
rollers themselves. The errors introduced by this problem can be removed by 
measuring the rolling slip at very low rolling speeds. 

2-2-4 THE MEASUREMENT OF THE TOROID SURFACE TEMPERATURE. 

In the analysis and reduction of the test data it is important that the disc 
temperature be used as the reference inlet temperature for the film thickness. This 
temperature can be measured by embedding a thermocouple directly below the surface of 
the toroid and then to take this signal out through mercury slip rings. This method 
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is however not very practical when a number of different toroids are involved and is 
also very costly from an installation point of view. With care the surface 
temperature can also be measured by using a trailing thermocouple that rides on the 
disc surface. The disadvantage of this technique is that it can be speed sensitive 
in its response because of frictional heating. 

The latter technique is employed here and care is taken to ensure that the contact 
force on the thermocouple is not excessive. A temperature reference bath junction 
ensures that the same reference level for the thermocouple is used at all times. The 
signal from the thermocouple is amplified using common mode rejection techniques to 
minimize the influence of electrical noise and other disturbances. Calibration is 
done by the boiling water method adjusted for sea level differences. This 
calibration is checked periodically. Only slight deviations are encountered. 

Because of the frictional heating at the junction/toroid interface a variation of 
about 2 °C is found in the signal between stationary discs and those rotating at a 
surface velocity of 120 m/s. The overall reproducibility of the temperature 
measurement is better than 2°C. 

The temperature on the machine is regulated through the use of heaters and 
coolers on the test fluid. This test fluid would be allowed to circulate freely 
before the start of a test series in order to bring the machine up to a uniform 
temperature. No specific effort is made however to maintain a given set point 
temperature during the test and. the reason for this will become clear in the analysis 
of the results. 

2-2-5 CONTACT RATIO MEASUREMENT. 

In order to measure the degree of asperity contact that took place at any one 
time an electrical conductance circuit as shown below is employed. 

The selection of the various electrical components is such that only very low 
voltages are applied across the two discs to minimize the risk of surface damage 
through arcing. The output signal of network is essentially DC in nature with only 
slow fluctuations in the contact conditions being recorded. 
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Schematic of the contact ratio measuring circuit. 

2-3 TRACTION MEASUREMENTS. 

Traction curves are obtained by the slow rotation of the upper assembly from a 
positive value of side slip/roll ratio to a negative value. The signals from the 
above discussed transducers are fed into a digitizer from where they are led into an 
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Apple 11+ computer for plotting and storage on magnetic media for future use. 

The computer would automatically trace the force versus slip curve on the screen. By 
reversing the direction of rotation of the machine, a duplicate set of curves can be 
obtained. For each experiment 500 data points are taken at fixed time periods of .2 
seconds. Multiple data points would be stored as separate entries. 

After the completion of a test series the data would be recalled into memory of the 
computer and further manipulated. This manipulation consisted of the averaging of 
the multiple entries, the filling in of any gaps in the data through forward and 
backward interpolation, the comparison of the traces for the forward and reverse 
rolling direction and the centering of the traces about the center lines. After the 
centering operation the data would be smoothed by a 'N 1 point averaging technique for 
traction points after the peak traction points. For storage a geometric series is 
used so that the total traction trace is now represented by 40 data points for each 
measured variable. These traces are then stored on magnetic media and used for 
further manipulation and data extraction at a later point 

2-3-1 TRACTION TEST RESULTS 

In order to shed some light on the observed efficiency trends as reported by 
Mitchell and Coy [1982] it was agreed that a select number of the fluids as used 
there would be tested for their traction. The selected fluids are; 


Designation in 
Mitchell & Coy 

F 

K 

E 

A 

I 


Commercial name 


Mobil RL-71 
Royco 555 
Syntech NS-6774 
Kendall ATF-220 
Mobil Jet II 


X Santotrac 50 (not in Mitchell and Coy) 

Y Sunoco TDF88 (not in Mitchell and Coy) 


Ideally these above lubricants should have come from the same batch as those used 
for the reported gear efficiency measurements. Unfortunately this was not possible. 

It should be stressed that while the designation may be the same, the actual fluid 
may not be. 

The viscosity temperature and the viscosity pressure behaviour of these fluids are 
shown in Figs 2-4 and 2-5. These results are from Present et al [1983], and will be 
used for the filmthickness calculations for the analysis of the various tests. The 
dependence of the fluid thermal conductivity on temperature is shown in Fig 2-6. 

While attempts were made to measure this property by an external source, it was found 
that the results showed too much uncertainty to be useful here. The dependency of 
the thermal conductivity on temperature was estimated by the method as used by 
Jackson [1980]. 
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To get adequate traction data on the above fluids for the purposes of traction 
prediction it was decided to test in under the following conditions; 


Mean Contact Pressure 
Mean rolling Velocity 
Mean inlet temperature 
Slip levels 


.7 - 1.3 GPa in 4 steps 
10, 20 and 30 m/sec 
30 and 80°C 
-.20 to .20 


Traction test on the first five fluids indicated above were taken under these 
conditions. The traction test results for the last two fluids were taken from 
earlier reported data. In total close to 200 traction tests are taken. These traces 
and a summary of all the test conditions may be found in Appendix I. 
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Fig. 2-6 


Thermal conductivity of the test fluids as a function of 
temperature. 
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3-0 THEORETICAL ANALYSIS 

In order to understand the required analysis of the experimental data it will be 
helpful to consider the following discussion of traction. 

The ability of a fluid film, trapped under high pressure in the elastically 
deformed region of two loaded curved elements, to transmit a tangential force froiri 
one element to the other, is commonly referred to as friction or traction. The 
magnitude of this force depends on several variables such as :1) the contact 
kinematic conditions of slip, spin and sideslip, 2) the fluid present, 3) 
temperature, pressure and operating speeds. We will examine the traction behaviour 
under simple slip only. 

Under conditions of increasing slip between the two elements an increasing 
traction force is transmitted up to a certain limit at which point it will decrease 
with further slip. See Fig. 3-1 



Fig. 3-1 Typical traction /slip curve. 

There are three regions identified on this traction curve and the behaviour 
in each of these regions can best be described by the Deborah number. For a simple 
Maxwell viscoelastic model this number is the ratio of the relaxation time and the 
mean transit time, see Johnson and Tevaarwerk [1977]. 

(A) The linear low slip region. Thought to be isothermal in nature, it is caused 
by the shearing of a linear viscous fluid (low De) or that of a linear elastic solid 
(high De). 

(B) The nonlinear region. Still isothermal in nature but now the viscous element 
responds nonlinearly. At low De this portion of the traction curve can be described 
by a suitable nonlinear viscous function alone, while at high De a linear elastic 
element interacts with the nonlinear viscous element. 

(C) At yet higher values of slip the traction decreases with increasing slip and 
it is no longer possible to ignore the dissipative shearing and the heat that it 
generates in the film. Johnson and Cameron [1967] showed that the shear plane 
hypothesis advanced by Smith [1965] does account for most of their experimental 
observations in this region. More recently Conry et al [1979] have shown that a 
nonlinear viscous element together with a simple thermal correction can also describe 

this region. 
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3-1 ISOTHERMAL TRACTION ANALYSIS. 

The rheological model that describes the traction under simple slip in all 
three regions of operation fairly well is the J & T traction model as presented by 
Johnson and Tevaarwerk [1977]; 

(3-D ?fr + F(T) = s 

G dt 

where T = shear stress [Pa] 

G = shear modulus [Pa] 

t = time [sec] 

g = shear strain rate [1/sec] 

The dissipative function F(T) is open to the choice of the researcher to fit the 
observed traction but Johnson and Tevaarwerk [1977] found that the hyperbolic 
sine; 

(3-2) F(T) = sinh 

v x s 

where T s = non-linear shear stress parameter. [Pa] 

V = local viscosity [Pa.sec] 

described all of their experimental results in regions (A) and (B) very well. At 
higher pressures and for fluids with high traction coefficients this dissipative 
function may be replaced by the purely plastic behaviour of the material; 

(3-3) F(T) = 0 for T < T c ; F(T) = g for T = T c 

where T c = limiting shear strength of the fluid [Pa]. 

Whether the perfectly plastic behaviour of the material is intrinsic is not clear. 

Work by Johnson and Greenwood [1980] suggests that it is possibly the result of 
thermal behaviour of the sinh model. For many applications the elastic/plastic 
traction model is adequate. It was used by Tevaarwerk and Johnson [1979c] and 
Tevaarwerk [1979b] to predict traction under various conditions of slip and spin. 

The analysis is completely isothermal in nature and for simple slip the traction is 
given by; 

(3-4) J4 = -§- ( tan-lS + (1 1 

2 J1 

where S = 

The shear strain rate in the fluid was taken to be the same everywhere in the 
contact and assumed to be constant throughout the film thickness. Its magnitude was 
taken to be; 



Equation (3-4) results from the integration of stresses, caused by the shearing 
of an elastic element of pressure independent average shear modulus G, and the 
plastic stresses proportional to the local Hertzian pressure. The predicted traction 
from an elastic/plastic model compares very well with the experimentally observed 
values for combinations of slip and spin, provided that the spin or slip are not too 
large. Large slip or spin results in almost purely dissipative stresses over the 
contact area and hence non-isothermal behaviour. Traction prediction under these 
conditions is still possible but the thermal effects need to be brought into the 
picture. Tevaarwerk [1981],[1979d] presents two techniques for calculating such 
spin traction curves. The latter technique requires the shape of the traction curve 
in the large slip regime to provide a simple correction to the isothermally predicted 
spin traction. 

Isothermal traction analysis was used in the reduction of traction data for the 
report on the results of the traction measurement, Tevaarwerk [1981b]. 

3-2 THERMAL TRACTION ANALYSIS. 

The ability to separate the elastic stresses from the plastic ones can now be used 
to perform a thermal traction calculation. The analysis presented here follows the 
technique outlined by Tevaarwerk [1980] and [1982]. 

Equation (3-4) resulted from the integration of isothermal elastic and plastic 
stresses over the contact area of an ellipse. For the region of contact under 
elastic stress the shear energy is conserved and therefore does not give rise to 
temperature increases. The plastically deforming region however, is non-conservative 
and a temperature rise in the fluid is expected. This will lead to a reduction in 
the local plastic strength of the fluid. Equation (3-4) may therefore be better 
written in its elastic and plastic parts; 


(3-5) 

J4 = J4p + J4 e 


(3-6) 

„ _ 4 S 

d4e ‘ ( 1+S 2 ) 2 


(3-7) 

2 

J4 n = - [ tan” 1 S + 
V t r 

S(S 2 -1) , 

(1+s 2 ) 2 J 


This modification is given by T c / T 0 where T c is the average 
stress under thermal conditions and T 0 is the average stress under isothermal 
conditions. We are dealing therefore with averaged stresses in the plastic region of 
the contact even though the isothermal stress distribution is according to the 
Hertzian pressure. This seemingly contradictory assumption is supported by 
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theoretical evidence by Tevaarwerk [1979d]. The modified equation would therefore 
be; 


(3-8) J4 t = J4 e + J4 p 

x 0 

The modification term T c /T 0 can be found from a thermal balance over the 
contact region under plastic stress. This region can be thought of as a thermal 
source whose heat is conducted/cohvected away. The length of the source is a 
function of the location of the onset of plastic deformation after the initial 
elastic region, however in this simple model we will take the source length to be 
"a" where this is the semi contact length in the running direction. As a simple 
thermal balance we will use the expression reported by Johnson and Cameron [1967] for 
the shear plane temperature; 


(3-9) 

®c “ ®o ~ *3 

h r « , l 

kf tFS+ VP e 1 

where; 

(3-10) 

P - 

r e j, 

K s 

1 - 

and 

q = average 

thermal strength of the source [W/m 2 ] 


0 C = shearplane temperature in the contact [°C] 

9 0 = inlet temperature of the fluid [°C] 

kf = thermal conductivity of fluid in the contact [W/m°C] 

h = central film thickness in the contact [m] 

k s = thenml conductivity of the roller material [W/m°C] 

C s = specific heat of the roller material [J/kg°C] 

p s = density of the roller material [kg/m 3 ] 

F s = thermal resistance of the film [-] 

This expression is only valid when the heat is conducted through the film and 
convected away by the discs, a condition that is true for most traction drive 
contacts. The strength q of the source is given by; 

(3-11) q = T AU 

The thermal resistance of the film F s may be calculated from the following 
expression by Tevaarwerk [1985]; 
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(3-12) F s = .07675 +.003455 Zj +.0000325 Zj 2 -.016 Ln(30/P e ) 

A + B P 
Where Zi = 8 0 + D 

In order to proceed any further we need a relationship between temperature and 
the shear strength of the fluid. A typical relationship that has its roots in the 
Eyring theory of fluid transport is given by; 

1 2 CE 

(3-13) T(0) = — [ A + BP + (0+D) Ln (- - )] 

where ; A - viscosity tenperature constant [°C] 

B = pressure viscosity constant [°C/Pa] 

C = non-linear shear stress constant [°C/Pa] 

D = fluid solidification tenperature [°C] 

E = fluid viscosity [Pa. sec] 

P = pressure of the fluid in the contact [Pa] 

At first sight it seems that this equation has five disposable constants in it, 
however two of these constants ( A and D ) may be obtained from the atmospheric 
viscosity temperature relationship, and one more ( B ) can be obtained from the Barus 
viscosity pressure relationship. The constant D is known as the solidification 
temperature; the temperature to which the fluid should be cooled to become solid like 
under atmospheric pressure. Only the constants C and E need to be determined 
experimentally from the traction results and this will be done in the next chapter. 

It should be noted here that the ultimate aim is to derive the fluid traction 
parameters such that they apply for all the experimental conditions reported here. 

By using equations (3-9), (3-10), (3-11) and (3-13) the average thermal shear 
stress can be obtained for a given set of conditions. In equation (3-8) we need the 
ratio, of the average contact shear stress under thermal conditions to that under 
isothermal conditions. This is really the ratio of the shear stress given by 
equation (3-13) evaluated at the shearplane temperature 6 C ( from equation 3-9) 
and the shear stress as evaluated at the inlet temperature conditions 9 0 . The 
other contact conditions remain the same for this ratio calculation. 

The method that is outlined above was used for the traction data analysis as 
reported by Tevaarwerk [1985] , with the exception of equation (3-12). Previously 
the value of F s was fixed at .1 . For the analysis of the traction data as 
reported here, several improvements were made in the model and these changes are 
outlined in the next chapter. 
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4-0 EXTRACTION OF THE TRACTION PARAMETERS. 

In order to use the traction model to calculate the level of traction in a given 
contact condition, the relevant traction parameters need to be extracted. In 
equation (3-1) we saw that both the elastic effects and the nonlinear viscous effect 
account for the behaviour of the fluid. Elastic effects are typically only felt 
under conditions of high contact pressure, low temperature and low slip. These 
conditions are not typical of those occurring in gear tooth contacts, so it would 
seem safe that we ignore any of the elastic fluid effects for the present analysis. 

This means that we only need to know the large strain viscous effects. 

4-1 EXTRACTION OF THE LARGE STRAIN PARAMETERS. 

The large slip region , that is the region beyond the traction peak, is 
exclusively governed by the dissipative element in the rheological equation. All the 
elastic effects have completely disappeared so that it is now quite easy to extract 
the governing parameters for this region. In essence what is required is a reverse 
analysis of the traction calculation normally used for the calculation of the 
traction curves. When the elastic response of the fluid is no longer dominant all 
the shear strain rate in the film is taken by the viscous element and we may express 
equation (3-1) as; 

(4-1) F(T) = I = -la- Sinh ( -1- ) 

v A s 

When the argument T/T s is larger than 1 this may be further simplified to the 
following; 

(4-2) F(T) = g = Exp ( — ) 

s 

This expression relates the stress to strain response for a small element of fluid 
that is subjected to a shear strain rate g , is at a given pressure and a given 
temperature. Each of the parameters in this equation can be a function of 
temperature or pressure or both. In the analysis of traction it is normally 
sufficient to treat the fluid viscosity as a function of temperature and pressure and 
to make the non linear stress parameter a function of temperature only. The following 
relationships will be used here: 

— for the viscosity we will combine the Vogels temperature viscosity and the Barus 
pressure viscosity relationships as follows; 


(4-3) 


V(0,p) = V Q exp ( 


A 

0+Dv 


+ 


pB 


) 


where V 0 = viscosity constant [Pa.sec] 

A = temperature viscosity constant [°C3 

D v = solidification constant for Vogels equation [°C] 
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B = pressure viscosity constant [°C/Pa] 

Dp = solidification constant for Barus equation [°C] 
6 = local temperature [°C] 
p = local pressure [Pa] 


— for the temperature dependence of T s on the local temperature we will use the 
following; 


C = non linear shear stress constant [°C/Pa] 

D s = solidification constant for Eyring equation [°C] 

These parameters ( except for the temperature and the pressure) are all constants 
for a given fluid. 

In order to proceed any further we have to make some assumptions about the 
conditions within the contact. Firstly we will assume that average conditions exist 
for the contact. Hence temperature and pressure are are only mean temperatures and 
pressures. Also we will assume that the important temperature for the analysis is the 
shear plane temperature as given by equation (3—9). Furthermore the shear strain 
( rate g will be assumed to be constant throughout the film and in the case of 

the traction experiments it is constant over the contact area. Its magnitude is 
given by; 

(4-5) g = Av/h [sec" 1 ] 

where h = central film thickness [m] 

The central filmthickness is calculated from the expressions by Hamrock and Dowson 
[1977]. Further modifications to this film thickness to allow for inlet shear 
heating are made by using the Murch and Wilson [1975] approach. 

Hence in total we will make the following simplifications; 


-6) 

: 

local 

average 

equation # 

temperature 

0 

®c 

(3-9) 

shear strain rate 

o 

% 

Av/h 

(3-4b) 

pressure 


p 

P=2P 0 /3 

(2-1) 

shear stress 


T 

T = F x /F z 



(4-4) 

where 
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By substituting the expressions from equation (4-3) to (4-6) into equation (4-2) we 


obtain; 

(4-7) 


Av _ 1 ( T C _ A _ B P 

h 2CVo(0 c +D s ) 0 C +D S 00+Dy 0^+Dp 


At a first glance there appear to be 7 parameters that can be used to fit the data. 
This is not the case however because several of these are determined by the results 
from other experiments. For example by the correlation of the viscosity data for the 
lubricants to equation (4-3) solves for the constants A, B, D v , V 0 and 
Dp. It remains therefore to find the constants D s and C. These can be 
derived from the thermal region of the traction curve itself by curve fitting 
equation (4-7) to it. For this purpose it is better to write this equation in a 
slightly different form; 


(4-8) 


T _ _1_ , A + BP + 2CVo Av 

e c +D s “ c e c +Dv e c +Dp (e c +D s) h 


From the above equation it is apparent that if this relationship holds then the 
results from the traction measurements should form a straight line when plotted in 
the above fashion. The slope of this line gives the value for C while the value of 
D s can be calculated by adjusting it until the best fit is obtained. 


4-2 REDUCED PRESSURE EFFECTS. 

When a fluid film is present in the contact zone the pressure distribution is not 
strictly Hertzian but modified by the fluid. Due to the hydrodynamic action the 
pressure distribution is more peaky, and is spread over a broader area. This will 
lead to a reduction in the mean pressure in the contact. The reduced pressure may be 
calculated from the pressure ratio as given below; 

(4-9) P P = P [1 - 4 G e -- 25 YI~ 3 e (‘ 2 * 3 / k ) 3 

where Gg = Johnsons elasticity constant [-1 
P = mean Hertz contact pressure [Pa] 

YI = Inlet shear heating factor [-3 
k = aspect ratio (b/a) [-3 

P r = reduced contact pressure [Pa] 


4-4 TRANSIT TIME DEPENDENCE. 

The fluid in the contact is subjected to very high pressure gradients for a very 
short time duration. It is therefore not likely that the viscosity in the contact 
has reached equilibrium conditions. How near equilibrium that it is will depend on 
several factors such as the time it is subjected to the pressure and the easy with 
which the molecules can rearrange themselves. 
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A suitable dimensionless grouping that takes these various factors into account 
given as; 


(4-10) 


* = JLXiM 

T a P 


is 


where § = dimensionless time delay parameter [-] 

U = mean roller speed [m/sec] 

V(0 o ) = viscosity at inlet temperature [Pa.sec] 
a = semi Hertz contact length [m] 

P 0 = Hertz contact pressure [Pa] 

This parameter is used in conjunction with the viscosity constant in the following 
way; 

(4-11) Vo = E0 $ E2 [Pa. sec] 


where E0 = equilibrium viscosity [Pa.sec] 
E2 = time delay exponent [-] 


The single constant V 0 has now been replaced by this two constants equation so 
we have increased our degree of freedom by one for the system. 

4-5 THERMAL CONDUCTIVITY TEMPERATURE DEPENDENCE 
As shown in Fig 2-6 the thermal conductivity of the fluid is a function of the 
temperature. This has a direct influence on the thermal resistance of the fluid in 
the contact. In order to take this effect into account the thermal conductivity was 
fitted to the following equation; 


(4-12) 


k f 


k’ 

e c + k" 


Since the shear plane temperature occurs in equation (4-12) and the thermal 
conductivity occurs in equation (3-9) the solution to both the thermal conductivity 
and the shear plane temperature must be done on an iterative basis. 


4-6 ANALYSIS OF THE EXPERIMENTAL RESULTS. 

By combining the above equations (4-8), (4-9) and (4-10) the traction data may be 
fitted to the following equation; 




0 c + Dv 


BP r . , 2 C EO^ 2 Av 

9 C +Dp L ( 8 C +D S ) h 


(4-13) 


T 

®C + ^S 
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where the shear plane temperature is calculated from equation (3-9) as; 


(3-9) e c - 9 0 = T Av [ Fs + 1 

and P e is given by; 

<«•> -• ■ ^ 

and F s is; 

(3-12) F s = .07675 +.003455 Zj +.0000325 Zj 2 -.016 Ln(30/P e ) 


the dimensionless temperature Zj is modified somewhat now to take the two 
different solidification temperatures into account; 


(4-11) 


_A_ + 

00+Dy 0o+D^)| 


From the experimental traction data as reported in Chapter 2 the curves where no 
asperity contact took place were used in a multiple curve regression analysis to 
determine the parameters C, E0 and E2. A summary of the traction curves used in the 
regression and the resulting fit of the traction data to equation (4-10) is shown in 
Appendix II for the 7 fluids investigated here. 

The regression constants for each of the fluids are given at the end of the summary 
sheet for each fluid. The variable El corresponds to C in the above equations. 
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5-0 PREDICTION OF THE GEAR MESH FRICTION LOSSES. 

The action of gear teeth that are transmitting power from one gear to the next may be 
modelled as a pair of rollers in contact. This modelling has allowed us to use 
traction test machines so that the essential rolling and slip traction data and 
models can be extended to the prediction of the gear mesh losses. What we have done 
so far is to take the traction data and to reduce all of it to a few fixed 
parameters through the traction modelling. With this model we should now be able to 
predict the level of traction for a given fluid when the contact pressure, rolling 
velocity, slip velocity and temperatures are known. While in the traction testers 
the parameters such as load, speed and contact radii are essentially constant during 
a given test, this is by no means the case in the contact of gear teeth, so we need 
to develop the proper equations for these first. After that we can then apply the 
traction prediction technique to a typical set of gears. 

In this chapter we will gradually work towards this goal and demonstrate how this 
is done for spur gears. The technique will only be demonstrated for spur gears here 
but is not limited to these alone, since it contains all the essential details needed 
for friction prediction in other types of gears. 

5-1 CONTACT CONDITIONS IN A PAIR OF SPUR GEARS 

In order to predict the local traction- at the contact formed between two gear teeth 
in load bearing we need to know the local curvatures and the velocities. These can 
be obtained from the basic gear parameters such as the operating centre distance, the 
gear ratio and the operating pressure angle, provided that the gear form is of the 
involute type. 

Consider a pair of spur gears as shown below. 
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If we assume that the tooth profile is of the involute form over the region that 
contact takes place than we may write the following relationships; 


(5-1) C d = ( Rb p + Rbg )/cos © 


where © = operating pressure angle [rad] 

C d = operating centre distance [m] 

Rbp = base radius of the pinion [m] 

Rbg = base radius of the gear [m] 


The instantaneous contact radii for the gear and pinion to the mesh point are given 
by; 


(5-2) instantaneous radius for gear = Rbg tan © - x 

instantaneous radius for pinion = Rbp tan © + x 

where the location of the mesh point is measured a distance x away from the pitch 
point in the direction of motion of the gear teeth. From these we may calculate the 
equivalent contact radius as. used in the Hertz contact calculations; 


(5-3) 


Cd sin © 
' (I+l ) 1 


[( I - X )( 1 + X )] 


where Rx e = equivalent contact radius in the rolling direction [m] 
I = Ng/Np = gear ratio [-3 

Ng = nunber of gear teeth [-] 

Np = nimber of pinion teeth [-] 

and X = — = dimensionless mesh point distance [-3 

C d sin © 


This value of the equivalent contact radius may now be used to calculate the Hertz 
contact pressure P 0 and the contact dimension f a f in the rolling direction. 

These quantities are needed in equations (4-9) to calculate the reduced contact 
pressure P P and for calculating the Peclet number in equation (3-10). Also needed 
for this is the normal load F£ on the tooth. This load depends on the location 
of x and other gear parameters. It will be assumed here that its magnitude has been 
established by the user. Assuming that the teeth have no crowning so that line 
contact will occur, we write the Hertz equation for contact pressure and size as 
follows; 
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(5-4) 

(5-5) 


and 


a = 


Po = 


VI 


5-2 INSTANTANEOUS ROLLING AND SLIDING VELOCITIES. 

The thickness of the film that is formed in between the two teeth in contact 
depends upon the local conditions of viscosity and surface velocity. Also the amount 
of sliding velocity at any given instant is needed to calculate the level of the mean 
shear strain rate in the film. Hence, needed are the instantaneous rolling velocity 
and the velocity difference at the contact location. Assuming that we have rigid 
body rotation for the two gears these velocities are given as; 


(5-6) U = w p Q [ 1 + X(I-1)/2I ] 

and 

(5-7) Up - U g = ( X w p C d sin © )/I 


where 


Up = rolling velocity for the pinion 
Ug = rolling velocity for the gear 
U = mean rolling velocity (Up+Ug)/2 
Wp = angular velocity of the pinion 


[m/sec] 

[m/sec] 

[m/sec] 

[rad/sec] 


5-3 GEAR TOOTH TEMPERATURES. 

In all the equations dealing with the shear plane temperature predictions we used the 
inlet temperature of the fluid. This is really the temperature of the bodies in 
contact in the immediate vicinity of the contact. In the traction tests this 
temperature was measured with the riding thermocouple so that the traction results 
were corrected for this on a point by point basis. In the case of gears however we 
normally have no direct control over this temperature, but rather control the fluid 
inlet temperature to the transmission itself. The tooth temperature, which 
corresponds more directly to our inlet temperature 9 0 , is determined by the 
amount of thermal resistance that it encounters in sheding the heat from the 
frictional losses at its surface. The exact determination of this resistance is too 
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complicated fop analysis here. Suffice it to say that a very good estimate 
of this inlet temperature may be obtained by the following; 

(5-8) 0 O = 0j + K t Q 

where 0i = inlet temperature of the fluid into the gearbox [°C3 

Kt = thermal resistance of gear tooth to ambient [°C.sec/N.m] 

Q = total amount of heat generated during a mesh cycle [N.m/sec] 

Since the amount of heat generated during a given mesh cycle depends on the traction 
of the fluid, which in turn depends on the local tooth temperature, it is easy to see 
that an iterative type solution is required to obtain 0 O . 

5-4 FRICTION EQUATION 

In chapter 4 we analyzed the traction data with the hyperbolic sine as the principle 
non-linear viscous element. In order to extract the fundamental constants from the 
traction data we had to make some assumptions about the magnitude of the term 
T/T s . When this term was much bigger than unity we approximated the hyperbolic 
sine with the exponential form. This allows us to use very simple regression 
techniques to find the fundamental constants for the given fluid by using equation 
(4-13). When dealing with traction drives where the traction is such that this 
stress ratio T/T s is always greater than 1 we can use equation (4-13) or the 
simpler form (3-13) directly to predict the traction as well. It is only when the 
slip velocities are high and when the traction properties of the fluid are such that 
we might violate the stress ratio condition , that we have to return to the original 
form of the Eyring non-linear viscous function. What this really means is that we 
have to substitute the traction constant that we solved for into equation (4-1). If 
we don’t do this we will not be able to predict the linear thermal viscous regime of 
the traction curve correctly. 

Starting now with equation (4-1) and backsubstituting the constants that we solved 
for we may write the following equation for the contact shear stress; 


(5-9) 


T = ( 


^c + Ds 

El 


) Sinh 


-i p_E 


B P T 


E0^2 exp ( OT 


0„ + D c 


- 


w 
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5-5 OUTLINE OF THE CALCULATION PROCEDURE. 

We are now in a position to calculate the traction (or friction) from the equations 
developed so far. Because of the nature of the process we will indicate the 
variables that are solved for in a form of flowchart. This flowchart will consist of 
three levels of calculations. These are; 

1) variables only dependent on the geometry 

2) variables dependent on the above and 0 O 

3) variables dependent on the above and 0 C 

Since there are two levels of iterations required at steps 2 and 3 these will 
become the inner most loops. Table 5-1 shows the flowchart that outlines the level 
at which each variable is calculated and the appropriate equation to be used. 
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Fixed parameters 


Gearing data : C^, ©, I, Wp > Fz» ^C> ^t 

Material data : E’, k s , C s , s 
Fluid data : A, B, VO, D v , Dp, K’, K", 0j 
Traction data : D s , EO, El, E2 


Variables dependent on geometry alone 


Symbol: 

Ryo 


Equation: 


Up-Ug 


Variables dependent on 0 O 


Symbol: 

Equation: 

®o 

5-8 

t 

4-10 

h 

Hamrock-Dowson 

Pr 

4-9 

0 

g 

4-5 

Zi 

4-11 

Pe 

3-10 

Fs 

3-12 


Variables dependent on 9 C 


Symbol: Equation: 
0- 3-9 


Table 5-1 : Flowchart of the traction calculation and the variables associated. 
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C" 5-6 EXAMPLE TRACTION PREDICTION. 

•As an example of the traction predicted by the method as outlined in Table 5-1, it 
was decided to do this for each of the fluids that were used in this investigation. 
Ultimately such a traction prediction would be an essential step in 
the prediction of the overall powerloss in a given transmission. Since we have not 
concerned ourselves with losses other than those caused by the sliding action of the 
gears it does only predict the sliding losses. Depending on the conditions in the 
transmission these losses may dominate or simply be a part of the overall. 

Here we will try to predict the powerloss for the fluids tested in traction and to 
see how these predicted losses compare with those measured by Mitchell and Coy 
[1982]. Since none of the gearing details are known to us we will have to make some 
assumptions about the geometry needed in our calculations. For simplicity we will 
make these predictions at the following conditions; 

R xe = .025 [m] 

U =10 [m/sec] 

(5-10) Po = 1 [GPa] 

Oj = 80 and 100 [°C] 

Up-Ug — > 0 to U [m/sec] 

Kt = 0 [°C.sec/N.m] 

v ... As a first prediction the amount of sliding traction is shown in Fig 5-1. This is for 

all the fluids tested here. The traction is plotted as the traction coefficient here 
which is defined as; 

(5-11) Traction coeff. = T/P = F x /F z 

These curves show the level of traction at a given value of slip. They do reveal that 
the traction for fluids X and Y are much higher than for the others. The next 
highest traction is revealed by fluid B. This order corresponds to the order in 
which the efficiencies measurements of Mitchell and Coy [1982] went. The curves in 
Fig 5-1 look like ordinary traction curves to us but this is misleading. Since only 
the viscous model is retained for the traction prediction done here, the initial 
slope is viscous in nature and not elastic as might be expected for some of the 
fluids. All this is not very important in the prediction of gear efficiencies 
because we want the total slip losses in the mesh. The instantaneous powerloss at a 
given value of slip is the product of the traction force and the slip velocity. In 
dimensionless form this is given as; 

(5-12) Instantaneous loss =%(~M » 

*z u p u g 



The curves of the instantaneous loss for the seven fluids are shown in figures 5-2 
and 5-3 for the two different temperatures. Again the order of the instantaneous 
losses are as those observed in the full transmission test particularly for the 
traction fluids at the higher temperature. These show a higher loss at the higher 
temperatures as was observed in the real transmission test. 

A yet more indicative parameter that shows the power loss in a mesh cycle would be 
the total of all the power loss from the starting slip till the current value of the 
slip. Essentially the area under the Instantaneous loss curves of figures 5-2 and 
5-3. These curves are shown in figures 5-4 and 5-5 for the two different inlet 
temperatures. 

5-7 THE CONSTANT FRICTION COEFFICIENT CONCEPT. 

Some of the gear friction models that are used in the prediction of mesh losses 
employ a constant friction coefficient. This friction coefficient provides for an 
easy method to calculate the losses if so desired. It also gives a very simple method 
for estimating the effect of the sliding friction loss on the overall losses of the 
powertrain. In order to use this method however we must have access to the average 
friction coefficient. This value may be calculated by the methods presented above 
because it is nothing more than the integral of the instantaneous losses divided by 
the current level of slip. For the example data as used above, these resulting 
constant friction coefficient curves are shown in figures 5-6 and 5-7. At 
sufficiently large values if slip it is indeed clear that the friction coefficient by 
this method is about constant. 
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Instantaneous Loss Fx/Fz#DU/U 
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Fig. 5-3 Dimensionless instantaneous loss factor for the test fluids as a 
function of slip. Met temperature = 100 °C. 
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***************** ******************************************'** 2 . 1 1360752E+ 1 1 
TRACTION CURVE DATA SUMMARY FOR M0BIL2 :Prog PLOTTER Rev# 1 
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Fluid- M0BIL2 : Po-1 GPa : Temp-31 'C 
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TRACTION CURUE DATA SUMMARY FOR NS-6774 ; Prog PLOTTER Rev# 1 
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ORIGINAL PASS II 
OF POOR QUALITY 

TRACTION CURVE DATA SUMMARY FOR ATF220 = Prog PLOTTER Rev# 1 
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Fluid' 5 RTF220 : Po=l GPa : Temp=26 'C 
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~* 10 r FI u i d*= RTF220 : Po«1.92 GPa : Temp*=3l 'C 
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~ ,10 n Fluid- RTF220 : Po-1.59 GPa : Temp~82 'C 
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TRACTION CURUE DATA SUMMARY FOR R0YC0555 = Prog PLOTTER Rev# i 
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10. 

7 

78 

200 

1 . 

7 

.00024 

R0YC0514ER 

1 

.01 

20. 

3 

73 

200 

1 . 

7 

.00024 

R0YCQ514FR 

1 

.01 

29. 

9 

79 

200 

1. 

.7 

.00024 

R0YC0514GR 

1 

.27 

10. 

3 

79 

400 

1. 

,7 

. 00030 

R0YC0514HR 

i 

.27 

20. 

3 

79 

400 

1 . 

7 

. 00030 

R0YC0514IR 

1 

.27 

30. 

. 1 

80 

400 

1. 

7 

.00030 

R0YC05 14 JR 

1 

.60 

10. 

,4 

79 

800 

1 . 

7 

.00037 

R0YC0514KR ' 

1 

.60 

20. 

, 3 

79 

800 

1. 

7 

. 00037 

R0YC0514LR 

1 

.60 

30. 

, i 

79 

300 

1 , 

7 

. 00037 

R0YC0514MR 

1 

. 93 

10. 

,6 

78 

1400 

1 , 

7 

.00045 

R0YC0514NR 

1 

.93 

20, 

n 

» 4 . 

79 

1400 

1 , 

7 

. 00045 

R0YCQ5140R 

1 

.93 

30, 

, 1 

81 

1400 

1 , 

7 

. 00045 



Traction Coeff Fx/Fz (-) Traction Coeff Fx/Fz C~) Traction Coeff Fx/Fz 
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/% 

l 


■ 1 0 r~ 

j 

.09 \- 

. 08 i- 

,07 b 
, 06 
,05 j— 
, 04 j- 
.03 j- 


Fluid* R0YC0555 


Po c l GPa : Temp*=31 'C 



.20 


10 P Fluid* R0YC0555 ; Po«1.2G GPa : Temp=35 'C 

. 09 \- 


.03 j— 


.04 f- 



• 10 [~ Fluid* R0YC0555 : Po*1.5S GPa : Temp=36 'C 

.0Sf- 

i 

.08 j- 

. 071 - 


.05 j- 




Traction Coeff Fx/Fz <-) Traction Coeff Fx/Fz ( ~) Traction Coeff Fx/Fz 
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. 10 


1 .09*- 

W * i 

.eat- 
.0?j- 

I 

.0Bp 

I 

, 05 H 


r Fluid- R0YC0555 : Po-1.92 GPa : Temp-37 'C 



r 


10 

03 j- 

08 h 

.0? f- 

,06F 
. 05 j— 
,04- 
, 03 j- 


Fluid- R0YC0555 : Po-1 GPa : Temp-78 'C 



, 10j — 

,09 i- 

i 

i 

,08r 


r 


07 

,06 

,05'F 

, 04 r 

t 

, 03 f* 


Fluid- R0YC0555 : Po-1. 26 GPa : Temp-79 'C 




Traction Coeff Fx/Fz (-) Traction Coeff Fx/Fz 
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Fluid" R0YC0555 : Po«1.59 GPa : Temp "7 8 'C 

i..09 r 

i 

.00p 

.07}- 

. 061 - 


.04 i- 



* I0 r Fluid" R0YC0555 ; Po«=1.92 GPa ! Temp"77 'C 

.09 h 

f 


.08 r 

i 

. 07 F 

i 

( 

.06 r 

i 

.05 1 — 


.04 



U*-10 m/'s 

U*“20 m/s 
U-30 m/s 


j. i i- 

.16 .18 
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TRACTION CURUE DATA SUMMARY FOR TDF53 =Prog PLOTTER Rev# 1 


Test # 

Po 

Uo 

To 

Fz 

kk 

a 

< ) 

( GPa ) 

( n / s ) 

( 5 C ) 

<N ) 

< - > 

( m ) 

TDF882B J8R 

i— ^ 

cn 

CD 

21.3 

57 

400 

.9 

.00036 

T0F882BJCR 

1.68 

30 . 0 

53 

400 

.9 

.00036 

TDF882BJDR 

1 .68 

40.5 

60 

400 

.9 

. 00036 

TDF882BJBR 

2.11 

22.0 

63 

800 

.9 

. 00046 

TDF882B J7R 

2.11 

30 . 0 

67 

300 

.9 

.00046 

TDF882BJ8R 

2.11 

40.5 

69 

800 

.9 

. 00046 

TDF382BJIR 

2.55 

21.8 

67 

1400 

.9 

. 00055 

TDF382BJ2R 

2 . 55 

30.4 

76 

1400 

.9 

.00055 

T0F882B J3R 

2.55 

40.5 

84 

1400 

.9 

.00055 

TDF832BJGR 

cn 

CD 

21.3 

79 

400 

.9 

.00036 

TDF982BJHR 

1 .68 

30 . 3 

32 

400 

a 

« W 

.00036 

TDF3S2BJIR 

1 .88 

40.5 

85 

400 

.9 

. 00036 

TDF882B JLR 

2.11 

23.4 

100 

300 

.9 

. 00046 

TDF882B JMR 

2.11 

31.3 

97 

300 

.9 

. 00046 

TDF882BJNR 

2.11 

40.5 

98 

800 

.9 

. 00046 

TDF882BJQR 

2.55 

22.3 

106 

1400 

.9 

.00055 

TDF382BJRR 

2.55 

29.7 

107 

1400 

.9 

.00055 

TDF3S2BJSR 

2.55 

40.5 

112 

1400 

.9 

.00055 




Traction Coeff Fx/Fz (-) Traction Coeff Fx/Fz (~) Traction Coeff Fx/Fz 
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Fluid" TDF88 s Po“1.67 GPa : Temp“57 'C 




U-21 m/s 

U-29 m/s 
U-40 m/'s 


.04 .86 .88 .10 .12 .14 .16 

Slide/Roll DU/U 


Fluid" TDF88 ; Po-2.11 GPa : Temp“68 'C 


U-22 m/'s 


U-29 m/'s 
U-40 m/s 


2 .04 .06 .08 .10 .12 .14 .16 

Slide/Roll DU/U C-> 

Fluid" TDF88 : Po=2.54 GPa ; Temp«67 'C 


U-21 m/'s 

U-30 m/'s 
U-40 m/s 


02 .04 .06 .08 .10 .12 .14 .16 

Slide/Roll DU/U (-) 



Traction Coeff Fx/Fz <~) Traction Coeff Fx/Fz <-) Traction Coeff Fx/Fz 
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SHde/Roll DU/U (-) 



SI 1 de/Ro 1 1 DU/U C-) 



SHde/Roll DU/U (-) 
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ft 


♦ a***********#****************-******************************* 2 . 1 lo5807o6E+l 1 

TRACTION CURUE DATA SUMMARY FOR SAN50 : Prog PLOTTER Rev# 1 
#************#**#^*******#*******■*#************************'**'******20 : 50 l .1 


Test # 

Po 

Uo 

To 

Fz 

kk 

a 

( ) 

( GPa ) 

< r/s ) 

( ’C ) 

( N ) 

(-) 

( R ) 

SAN5044K1R 

2 . 55 

11.9 

37 

1400 

.9 

.00055 

SAN5044K2R 

2.55 

22 . 9 

50 

1400 

.9 

.00055 

SAN5044K3R 

2.55 

33 . 3 

57 

1400 

.9 

.00055 

SAN5044K4R 

2.11 

11.8 

38 

800 

.9 

. 00046 

SAN5044K5R 

2.11 

23.4 

40 

800 

.9 

. 00046 

SAN5044KGR 

2.11 

32.5 

44 

800 

. 9 

. 00046 

SAN5044K7R 

1 .S3 

32.7 

41 

400 

. 9 

. 00036 

SAN5044K8R 

07 

CO 

21.4 

37 

400 

.9 

.00036 

SAN5044K9R 

1 . 68 

11.4 

34 

400 

.9 

. 00036 

SAN5044JDR 

1.68 

11.8 

69 

400 

.9 

. 00036 

SAN5044JER 

1.68 

21.3 

67 

400 

.9 

. 00036 

SAN5044 JFR 

1.68 

31 . 1 

69 

400 

.9 

. 00036 

SAN5044JGR 

2.11 

12 2 

•y n 

800 

.9 

. 00046 

SAN5044 JHR 

2.11 

22.8 

74 

800 

.9 

. 00046 

SAN5044 J IR 

2.11 

34.0 

73 

800 

.9 

. 00046 

SAN5044J JR 

2.55 

32.7 

87 

1400 

.9 

.00055 

SAN5044JKR 

2.55 

22.9 

90 

1400 

.9 

.00055 

SAN5044JLR 

2.55 

11.4 

91 

1400 

.9 

.00055 
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Slide/Roll DU/U C-> 



Slide/'Roll DLJ/U (-) 



Slide/'Roll DU/U (-) 


. 18 


_i 

. 18 


. 18 



Traction Coeff Fx/Fz (-) Traction Coef-F Fx/Fz ( ~ ) Traction Coeff Fx/Fz 
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Slide/RoU D’J/U <-) 



.06 .08 .10 .12 
Slide/Roll DU/U (-) 



,Ub . uo . i kj .»c 

SlIde/RoU DU^U (-) 







Traction Coeff Fx/Fz (-) Traction Coeff Fx/Fz (-) Traction Coeff Fx/Fz 
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Traction Coeff Fx/Fz (-) Traction Coeff Fx/Fz 
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************************************************************* 2.11358073EE+11 
MULTIPLE TRACTION CURUE DATA ANALYSIS FOR M0BIL2 ; Prog HYPER Rev# 7.4 


Data corrections: 

Long, slip YES: Reduced press. YES 


Inlet heating YES: Therm. Cond.YES 


The following oil properties were used in the analysis: 


Uiscosity temperature 
Uiscosity Pressure 
Thermal Conductivity 
Eyring Temperatures 


Av= 378.0 ’C Dv=120 . 0 ’C 00= .0000853 Pas 
B=2448 . 0 ’C/GPa Dp=174.0 ’C 

K ’ =1851 . 0 W/m’C K ’ ’ =272 . 0 ’C 
Dp 174.0 ’C Ds =120.0 ’C 




TEST INCLUDED IN THE MULTIPLE TRACTION DATA ANALYSIS 


Test # 

Po 

Uo 

To 

Fz 

kk 

a 

Pes 

( — ) 

( GPa ) 

( m / s ) 

( ’C ) 

(N ) 

(- ) 

( m ) 

' (- ) 

M0BIL5132R 

1 .01 

19.9 

33 

200 

1.7 

.00024 

5.31E+02 

M0BIL5133R 

1.01 

30.2 

35 

200 

1.7 

.00024 

8.84E+02 

M0BIL5134R 

1.27 

10.2 

34 

400 

1.7 

. 00030 

1 . 17E+02 

M0BIL5135R 

1 . 27 

20.3 

35 

400 

1.7 

. 00030 

3.73E+02 

MOB IL513BR 

1.27 

30.2 

37 

400 

1 .7 

, 00030 

6.23E+02 

M0BIL5137R 

1 .80 

10.3 

35 

800 

1.7 

. 00037 

7.87E+01 

M0BIL5138R 

1.80 

20.2 

38 

800 

1.7 

. 00037 

2 . 59E+02 

M0BIL5139R 

1.80 

29.1 

38 

800 

1.7 

.00037 

4.05E+02 

M0BIL513AR 

1.93 

10.9 

37 

1400 

1.7 

.00045 

6 . 06E+0 1 

M0BIL513BR 

1 . 93 

20.6 

38 

1400 

1.7 

.00045 

1 .8SE+02 

M0BIL513CR 

1 . 93 

27.8 

41 

1400 

1.7 

. 00045 

2.S8E+02 

M0BIL5149R 

1.80 

30.1 

83 

800 

1.7 

. 00037 

1 .31E+02 

M0BIL514BR 

1 . 93 

20.2 

84 

1400 

1.7 

.00045 

4 . 48E+01 


TYPE OF FIT IS 4 

BEST FIT CONSTANTS ARE: 

E0= 2.344SE-05 El= 1 . 9894E-05 E2= 3.5903E-01 E3= 0.0000E+00 
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c ( 0 d xB ^n x «^ 3 )Don x z 3 + <AQ+o)/ A y + (da+co/axd + ( (sq+oxqh )/na x 13x03x2)301 



4 

j 

-i 

-j 02 

1 

i 

j 

\ 

•4 

i 

1 

y y 23I80W JOJ- uo^oBvi^ j-o %\± eu^s o^oqjedXH 



Tau x El/CO+Ds) 
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#*********+#*##*********#*+********************************** 2 . 1 13580736E+ 1 1 
MULTIPLE TRACTION CURUE DATA ANALYSIS FOR NS-6774 ; Prog HYPER Rev# 7.4 


Data corrections; 

Long, slip YES : Reduced press. YES; 


Inlet heating YES = Therm. Cond.YES 


The following oil properties were used in the analysis: 


Viscosity temperature 
Viscosity Pressure 
Thermal Conductivity 
Eyring Temperatures 


Av= 793.0 ’C Dv=100 . 0 ’C 00= 
B=2788 . 0 ’ C/GPa Dp=140.0 ’C 
K ’ =1851 . 0 W/m’C K ’ ’=233.0 ’C 
Op 140.0 ’C Ds =100.0 ’C 


.0000982 Pas 




TEST INCLUDED IN THE MULTIPLE TRACTION DATA ANALYSIS 


Test # 

Po 

Uo 

To 

F 2 

kk 

a 

Pes 

( ) 

( GPa ) 

( m/s ) 

( ’C ) 

( N ) 

(- ) 

( m ) 

( - ) 

NS_675191R 

1.01 

10.2 

34 

200 

1.7 

.00024 

3.96E+02 

NS_G75193R 

1.01 

29.7 

38 

200 

1.7 

.00024 

1.70E+03 

NS_B75195R 

1.27 

20.3 

37 

400 

1.7 

. 00030 

7.81E+02 

NS_G75 1 96R 

1.27 

30.2 

40 

400 

1.7 

. 00030 

1 . 19E+03 

NS_675 1 97R 

1 .60 

10.2 

38 

800 

1.7 

. 00037 

1.84E+02 

NS_675 1 98R 

1 . S0 

20.3 

39 

800 

1.7 

. 00037 

5.33E+02 

NS_G75 1 99R 

1 .60 

30.2 

41 

800 

1.7 

. 00037 

8.44E+02 

NS_G75 1 9AR 

1.93 

11.1 

39 

1400 

1.7 

.00045 

1.4GE+02 

NS_S75i 9BR 

1.93 

21 . 1 

41 

1400 

1.7 

.00045 

3.97E+02 

NS_675 1 9CR 

1.93 

31.3 

45 

1400 

1.7 

.00045 

G.20E+02 

NS_675 1 SLR 

1 .60 

30.2 

81 

800 

1.7 

. 00037 

2.8GE+02 

NS_675 1 9NR 

1.93 

21.0 

84 

1400 

1.7 

.00045 

9.99E+01 

NS G75190R 

1 . 93 

30.5 

83 

1400 

1.7 

. 00045 

2.07E+02 


TYPE OF FIT IS 4 

BEST FIT CONSTANTS ARE: 

E0= 3 . 1248E-06 El= 2.9719E-05 E2= 1.7973E-01 


E3= 0 . 0000E+00 
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< (0d** )/'n x **3)D01 x 23 + (Aa+0)/Ay + ( da+0)/S x d + C ( sa+O ) X0H ) /HQ* l 3*03*2 ) D01 






t>X<19-SN 


joj. 


UO J. 


j-o tuj. aujs o^ioqjsdXH 



Tau x El/CO+Ds) 
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Of POOR 


**#****#****************##***+**************+*******+******** 2 . 1 1358073GE+1 1 
MULTIPLE TRACTION CURVE DATA ANALYSIS FOR ATF220 = Prog HYPER Rev# 7.4 


Data corrections: 

Long, slip YES= Reduced press. YES = Inlet heating YES: Therm. Cond.YES 


The following oil prop 
Viscosity temperature: 
Viscosity Pressure : 
Thermal Conductivity : 
Eyring Temperatures ; 


erties were used in the analysis: 
Av= 993.0 ’C Dv=120 . 0 ’C V0= 
B=3095 . 0 ! C/GPa Dp=182.0 ’C 
K ’ =1851 . 0 W/m’C K ’ ’ =252 . 0 ’C 
Dp 182.0 ’C Ds =200.0 ’C 




.0000650 Pas 


*********** 18 : 11 : 


52 


TEST INCLUDED IN THE MULTIPLE TRACTION DATA ANALYSIS 


Test # 

Po 

Uo 

To 

Fz 

kk 

a 

Pes 

( __ — ) 

( GPa ) 

( m/ s ) 

( ’C ) 

( N ) 

(- ) 

( m ) 

(- ) 

ATF224CV1R 

1.01 

10.2 

27 

200 

1.7 

.00024 

4 . 79E+02 

ATF224CV2R 

1.01 

20.2 

28 

200 

1.7 

.00024 

1.26E+03 

ATF224CV3R 

1.01 

30.2 

30 

200 

1.7 

.00024 

1.89E+03 

ATF224CV4R 

1.01 

20.2 

30 

200 

1.7 

.00024 

1 . 1 9E+03 

ATF224CV5R 

1.01 

10.0 

29 

200 

1.7 

.00024 

4.24E+02 

ATF224CV6R 

1.27 

9.8 

29 

400 

1.7 

. 00030 

3.19E+02 

ATF224CV7R 

1.27 

20.0 

30 

400 

1.7 

.00030 

8.73E+02 

ATF224CV8R 

1.27 

30 . 0 

31 

400 

1.7 

. 00030 

1 . 38E+03 

ATF224CV9R 

1.60 

10.2 

30 

800 

1.7 

. 00037 

2.15E+02 

ATF224CVAR 

1 .60 

20.2 

31 

800 

1.7 

. 00037 

6.36E+02 

ATF224CVBR 

1.60 

30.2 

33 

800 

1.7 

. 00037 

9.80E+02 

ATF224CVCR 

1.93 

10.9 

32 

1400 

1.7 

. 00045 

1 .42E+02 

ATF224CVDR 

1.93 

20.6 

33 

1400 

1.7 

.00045 

4 . 53E+02 

ATF224CVER 

1.93 

30.5 

35 

1400 

1.7 

. 00045 

7.37E+02 

ATF225 133R 

1 .01 

25.3 

84 

200 

1.7 

.00024 

3 .34E+02 

ATF225 134R 

1.01 

30.1 

85 

200 

1.7 

.00024 

4 . 41E+02 

ATF225 137R 

1.27 

30 . 0 

82 

400 

1.7 

. 00030 

3.44E+02 

ATF225 138R 

1.60 

10.5 

82 

800 

1.7 

. 00037 

3 . 06E+0 1 

ATF225 139R 

1.60 

20.3 

32 

800 

1.7 

. 00037 

1 .25E+02 

ATF225 13AR 

1.60 

30.1 

84 

800 

1.7 

. 00037 

2.S4E+02 

ATF225 13BR 

1 . 93 

20.9 

83 

1400 

1.7 

.00045 

9.40E+01 

ATF225 1 3CR 

1.93 

25.7 

83 

1400 

1.7 

.00045 

1 .40E+02 

ATF22513DR 

1 . 93 

10.9 

81 

1400 

1.7 

.00045 

2.58E+01 

TYPE OF FIT 

IS 

4 






BEST FIT CONSTANTS 

ARE: 






E0= 2 . 940 IE 

-06 EL 

= 3 . 343 

2E-0S 

E2 = 

4. 3173E-01 

E3= 0 . 0000E+00 



T au x El/(0+Ds) 


30 


r 


Hyperbolic Sine fit of traction data for RTF220 



LOG(2xE0xElxDU/CH0x(O+Ds)) + PxBXCO+Dp) + Rv/CO+Dv) + E2xL0G(EtaxU/C axP0) ) 
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# #***********#********************************************'*'** 2. 11058070 6E+11 
MULTIPLE TRACTION CURVE DATA ANALYSIS FOR R0YC05E5 ; Prog HYPER Rev# 7.4 


Data corrections: 

Long, slip YES: Reduced press. YES: Inlet heating YES: Them. Cond.YES 


The following oil properties were used in the analysis 


Viscosity temperature 
Viscosity Pressure 
Thermal Conductivity 
Eyring Temperatures 


Av= 8S5.0 ’C Dv=l 10.0 ’C V0= 
B=2823 . 0 ’ C/GPa Dp=225.0 ’C 
K ’=1851 .0 W/m’C K’ ’=221 .0 ’C 
Dp 225.0 ’C Os =110.0 ’C 


0000865 Pas 






TEST INCLUDED IN THE MULTIPLE TRACTION DATA ANALYSIS 


Test # 

Po 

Uo 

To 

Fz 

kk 

a 

Pes 

( ) 

(GPa ) 

( m/s ) 

( ’C ) 

(N ) 

( - ) 

( m ) 

( - > 

RQYC05142R 

1.01 

20.3 

34 

200 

1.7 

.00024 

3.42E+02 

R0YC05143R 

1.01 

30.1 

36 

200 

1.7 

.00024 

1 .26E+03 

R0YC05145R 

1.27 

20.3 

36 

400 

1.7 

. 00030 

5.77E+02 

R0YC0514GR 

1.27 

30.6 

39 

400 

1.7 

. 00030 

8.34E+02 

R0YC05147R 

1.60 

10.2 

36 

800 

1.7 

. 00037 

i .41E+02 

R0YC05148R 

1.60 

20.1 

38 

800 

1.7 

. 00037 

3.82E+02 

R0YC05 149R 

1.60 

30.1 

40 

800 

1,7 

. 00037 

6. 19E+02 

R0YC0514AR 

1.93 

10.5 

37 

1400 

1.7 

. 00045 

1 .04E+02 

R0YC05 14BR 

1.93 

20.4 

40 

1400 

1.7 

. 00045 

2.33E+02 

R0YC05 14CR 

1.93 

30.2 

43 

1400 

1.7 

.00045 

4.41E+02 

R0YC05 1 4LR 

1.60 

30.1 

79 

800 

1.7 

. 00037 

2 . 36E+02 

R0YC05 1 4NR 

1 . 93 

20.2 

79 

1400 

1.7 

.00045 

8. 19E+01 

R0YC05 140R 

1 . 93 

30.1 

81 

1400 

1.7 

.00045 

1 .66E+02 

TYPE OF FIT 

IS 

4 






BEST FIT CONSTANTS 

ARE: 






E0= 1.0050E- 

-06 El 

= 1 . 7647E-05 

E2 = 

4.982 

9E-01 

E3= 0.0000E+00 


ORICai^JAi. 

OF POOR QUALITY 
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( C0dx^ )/'fl x ^^3)D0"l x B3 + ( A d-f-0)/ A y + C^Q+0)/H x d + ( C s Q+0) x 0H)/ntf x I3 x 03 x 2)DCn 



1 




y 

SSQOOAOd JO J* uoiiobj* J-O auis o^oqjadAn 


J 

J 0E 



Tau x E 1 / ( O+Ds ) 
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MULTIPLE TRACTION CURUE DATA ANALYSIS FOR SAN50 : Prog HYPER Rev# 7.4 
Data corrections: 

Long, slip YES = Reduced press. YES ; Inlet heating YES : Them. Cond.YEs 

The following oil properties were used in the analysis: 

Uiscosity temperature: Av= S9S.0 ’C Dv= 65.0 S C U0- .000181G Pas 

Uiscosity Pressure : B = 3298.0 ’C/GPa Dp- G5.0 ’ .C 

Thermal Conductivity = K ,= 1506.0 W/m’C K ,,= 160.0 ’C 

Eyring Temperatures : Dp S5.0 ’C Ds - G5.0 ’C 


TEST INCLUDED IN THE MULTIPLE TRACTION DATA ANALYSIS 


Test # 

Po 

Uo 

To 

Fs 

kk 

a 

Pes 

( ) 

( GPa ) 

( m/s ) 

( ’C ) 

(N ) 

( - ) 

( m ) 

(- ) 

SAN5044 JDR 

1 .68 

11.8 

69 

400 

.9 

. 00036 

1 .94E+02 

SAN5044 JER 

1.68 

21.3 

67 

400 

.9 

.00036 

7. 12E+02 

SAN5044 JFR 

1.68 

31 . 1 

69 

400 

.9 

. 00036 

1.31 E+0o 

SAN5044 JGR 

2.11 

j_ 9 9 

72 

800 

.9 

.00046 

1 . 1 3E+02 

SAN5044JHR 

2 .11 

22.8 

74 

800 

.9 

. 00046 

4.52E+02 

SAN5044JIR 

2.11 

34.0 

78 

800 

.9 

. 00046 

9.42E+02 

SAN5044 J JR 

2.55 

32.7 

87 

1400 

.9 

'. 00055 

4.78E+02 

SAN5044 JKR 

2.55 

22 . 9 

90 

1400 

.9 

.00055 

1 .70E+02 

SAN5044 JLR ' 

2.55 

11.4 

91 

1400 • 

.9 

. 00055 

3.30E+01 

SAN5044K 1R 

2.55 

11.9 

37 

1400 

.9 

.00055 

4.97E+02 

SAN5044K2R 

2.55 

22 . 9 

50 

1400 

.9 

.00055 

9.21E+02 

SAN5044K3R 

2.55 

33 . 3 

57 

1400 

.9 

.00055 

1 . 33E+03 

SAN5044K4R 

2.11 

11.8 

38 

800 

. 9 

.00046 

6.81E+02 

SANS044K5R 

2.11 

23.4 

40 

800 

.9 

.00046 

1 .G8E+03 

SAN5044KGR 

2.11 

32.5 

44 

800 

.9 

. 00046 

2.45E+03 

SAN5044K7R 

1 .68 

32.7 

41 

400 

.9 

. 00036 

3.22E+03 

SAN5044K8R 

1.63 

21.4 

37 

400 

.9 

. 000o6 

2 . 1SE+03 

SAN5044K9R 

1.68 

11.4 

J 4 

400 

.9 

. 00036 

1 . 1SE+03 

TYPE OF FIT 

IS 

4 






BEST FIT CONSTANTS 

ARE: 






E0= 1.01G0E 

-03 El 

= 2.853 

3E-05 

i E2= 

2 . 7131E-0 1 

E3= 0 . 0000E+00 
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c ( 0 dx«)/n x ^* 3 )Dcn *23 + caq+o/ab + (da+o)/a><d + c c s a+o) *01-0 /na* 13x03*21301 
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************************************************************* 2.1 1 358 0736 E+l 1 
MULTIPLE TRACTION CURVE DATA ANALYSIS FOR TDF88 =Prog HYPER Rev# 7.4 


Data corrections: 

Long, slip YES ; Reduced press. YES: Inlet heating YES: Therm, Cond.YES 


The following oil properties were used in the analysis: 

Viscosity temperature: Av— S99.3 ’C Dv= 35.0 ’C V0- .0001489 Pas 

Viscosity Pressure : 8 5= 2975.0 'C/GPa Dp= 84. 0 ’C 

Thermal Conductivity : K’=1506.0 W/m’C K’ ,= 160.0 ’C 

Eyring Temperatures : Dp 84.0 Ds = 85.0 } C 


TEST INCLUDED IN THE MULTIPLE TRACTION DATA ANALYSIS 


Test # 

Po 

Uo 

To 

Fz 

kk 

a 

Pes 

( ) 

( GPa ) 

( m / s ) 

( ’C 

) (N ) 

(-> 

( m ) 

(- ) 

T0F8S2B J1R 

2.55 

21.8 

67 

1400 

.9 

.00055 

2.61E+02 

TDF882BJ2R 

2.55 

30 . 4 

76 

1400 

.9 

.00055 

3.84E+02 

TDF382BJ3R 

2.55 

40.5 

84 

1400 

.9 

.00055 

4.94E+02 

TDF882BJ6R 

2.11 

22.0 

69 

800 

.9 

. 00046 

3 . 38E+02 

TDF882B J7R 

2.11 

30.0 

67 

800 

.9 

. 00046 

6.24E+02 

TDF882E J8R 

2.11 

40.5 

69 

800 

.9 

.00046 

9.94E+02 

TDFS82B JBR 

1.68 

21.3 

57 

400 

. 9 

. 00036 

6.91E+02 

TDF882B JCR 

1 .S3 

30.0 

5-8 

400 

.9 

. 00036 

1.17E+03 

TDF882B JDR 

1.68 

40.5 

60 

400 

.9 

. 00036 

1 . 60E+03 

TDF882B JGR 

1.68 

21.3 

79 

400 

.9 

. 00036 

3. 19E+02 

TDF882BJHR 

1.68 

30 . 3 

82 

400 

. 9 

.00036 

5.58E+02 

TDF382BJ IR 

1 .68 

40.5 

85 

400 

.9 

. 00036 

8.67E+02 

TDF382S JLR 

2.11 

23.4 

100 

800 

.9 

. 00046 

1 .97E+02 

TDF8S2B JMR 

2.11 

31.3 

97 

800 

.9 

. 00046 

2.98E+02 

TDF882B JNR 

2.11 

40.5 

98 

800 

.9 

. 00046 

4 . SSE+02 

TDF882B JQR 

2.55 

22.8 

106 

1400 

.9 

.00055 

9 . 33E+01 

TDF882B JRR 

2.55 

29.7 

107 

1400 

.9 

.00055 

1 .52E+02 

TDF382B JSR 

2.55 

40.5 

112 

1400 

.9 

.00055 

2.79E+02 

TYPE OF FIT 

IS 

4 






BEST FIT CONSTANTS 

ARE: 






E0= 2.0053E- 

-03 El 

= 2.52 

49E- 

05 £2=- 

-1.770 

2E-01 

E3— 0.0000E+00 
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( C0dX B )/Tlx«*3)DCnx23 + (Aa + 0)/Ay + 

0E 03 

j — t j - i - t i • i t — t i — ■ I i r r • 


(da+e)/a><d + c csa+e)X0H i/na* 13x03x2)^01 

01 0 

I t »- t I { I I » I I 1 1 r \-~ja 
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*************************+*******************•**************** 2.1135807o6E+li 
MULTIPLE TRACTION CURVE DATA ANALYSIS FOR RL-714 = Prog HYPER Rev# 7.4 



Data corrections: 

Long, slip YES: Reduced press. YES: Inlet heating YES = Therm. Cond.YES 


The following oil properties were used in the analysis 1 


Viscosity temperature 
Viscosity Pressure 
Thermal Conductivity 
Eyring Temperatures 


Av= 75S.0 'C Dv=100 . 0 ’C V0= 
B=3363 . 0 ’ C/GPa Dp=210.0 'C 
K ’ =1851 . 0 U/m’C K ’ ’ =263 . 0 ’C 

Dp 210.0 ’C Ds =100.0 ’ C 


.0000367 Pas 




TEST INCLUDED IN THE MULTIPLE TRACTION DATA ANALYSIS 


Test # 

Po 

Uo 

To 

Fz 

kk 

a 

Pes 

( ) 

(GPa ) 

( m / s ) 

( ’C ) 

( N ) 

<-> 

( m ) 

<- ) 

RL_71S182R 

1.01 

20.3 

29 

200 

1.7 

.00024 

7.64E+02 

RL_715183R 

1.01 

30.2 

3 1 

200 

1.7 

.00024 

1 .20E+03 

RL_715184R 

1.27 

10.2 

30 

400 

1.7 

. 00030 

1 .82E+02 

RL_7 15 185R 

1.27 

20.2 

32 

400 

1.7 

. 00030 

5. 10E+02 

RL_7 15186 R 

1 .27 

30.2 

o3 

400 

1.7 

. 00030 

8. 18E+02 

RL__715 187R 

1.60 

10.2 

30 

800 

1.7 

. 00037 

1 .31E+02 

RL_715188R 

1.60* 

20.2 

32 

800 

1.7 

. 00037 

3.68E+02 

RL_715 189R 

1.60 

30 . 0 

35 

800 

1.7 

. 00037 

5.36E+02 

RLJ71518AR 

1.93 

10.6 

32 

1400 

1.7 

. 00045 

9 . 30E+0 1 

RL_71518BR 

1.93 

20.7 

35 

1400 

1.7 

.00045 

2.76E+02 

RL_715 18CR 

1.93 

30.5 

38 

1400 

1.7 

.00045 

4. 18E+02 

RL_7 15 1 SNR 

1 . 93 

20.9 

34 

1400 

1.7 

.00045 

5.48E+01 

RL 715130R 

1.93 

31 . 1 

84 

1400 

1.7 

.00045 

1 . 18E+02 


TYPE OF FIT IS 4 


BEST FIT CONSTANTS ARE: 

E0= 2.7050E-07 El= 2.4535E-05 E2= 4.9813E-01 E3= 0.0000E+00 




Tau x El/CO+Ds) 


30 


20 


Hyperbolic Bine fit of traction data for RL—714 




X 


x 


x 





x 




20 30 

LOG ( 2xE0xE 1 xDU/ C H0x ( O+Ds ) ) + PxB/CO+Dp) + Rv/CO+Dv) 4- E2xL0G(Et axll/{ axP0) ) 
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